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ABSTRACT 

A theoretical solution of the title problem was obtained 
and a computer program was developed to examine various geometric 
arrangements of elliptical holes in flat plates with different 
loading c o nulls i ons • 

Graphs of theoretical predictions are presented, discussed 
and compared with experimental results of other authors. 
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NOMENCLATURE 


2 C 


l 

Distance between the foci of the ellipses 


ct,b 


Major and minor semiaxes of the elliptical hole 


°7 


Elliptical coordinates of a point P(x,y) 


ro 


Constant value of f of any point on the elliptical 
hole 


$,Xi 


Stress functions 




Normal stress intensity in the tangential 
direction to =. constant 


<55 


Normal stress intensity in the normal direction 
to % = constant 


r $«f 


Shear stress intensity 


6* , By , Xxj 


The usual notation for stress intensities in 
Cartesian coordinates 


w* 


Dimension of the plate in the x direction 


Wy 


Dimension of the plate in the y direction 


}* = 

Vx 

q y = 

^ Wy 

6 l 


Uniformally distributed load applied on the 
edge x= and parallel to X axis 


Bz 


Uniformally distributed load applied on the 


a , b 


edge y= and parallel to y axis 

Semiaxes of the confocal ellipse ^-§0 passing 
through a point P’Cx^y) = P ( ^ > y 


SCF 


Stress concentration factor =. ( ^ rr,<xx 0 r 6rrr ' Co< 

6 *. 


Sc f' 


Stress concentration factor = SCF« C1 _ V) or 
SC F * Ci-},) 



- 6 - 



jj.= b/a. 



when loading parallel bo the minor axis of the 
elliptical hole. 



fJL - a/b 



when loading parallel to the major axis of the 
elliptical hole. 
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INTRODUCTION 



The behavior of thin flat elastic plates perforated with 
circular holes is -well known for in-plane loading. 

Theoretical methods and experimental results have been 
published by many authors. Refs. ( 1 - 7 ) 

The case of one elliptical hole in an infinitely large 
plate is discussed in reference [l] for several loading conditions 
including the case where the major axis of the ellipse is inclined 
to the axis of the axial load. 

Durelli and Murray [10] performed photcelastic tests on thin 
flat plates loaded with various combinations of biaxial stresses. 

Durelli, Parks and Feng [ 9 ] examined cases of rectangular 
plates with one elliptical hole under axial loading, again using 
the photoelastic method. 

The purpose of this study is to approach the problem 
theoretically and obtain results for the case of finite thin 
rectangular plates perforated with one elliptical hole and 
subjected to uniaxial or biaxial loading. 

A computer program was developed and the published results 
of available experiments were used for comparison. 



- 8 - 



THEORY 



The use of elliptical coordinates simplifies the mathematical 
analysis of the stresses around an elliptical hole in a flat plate 
loaded in the plane of the plate. 

The equations of transformation from Cartesian to elliptical 
coordinates are 



by a family of confocal ellipses and a family of confocal hyperbolae 
respectively, that are orthogonal and may be expressed as 



X = c cosK t COS 



( 1 ) 



y = csinh| silt ^ 



( 2 ) 



The conditions $ = constant and 'n = constant are represented 



c ' 1 ccsh^l 




i 



( 3 ) 





( J 0 



Figure la represents the elliptical coordinates ^ , 
of a point P <x,y ) . 
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Using the notation 



f~ JL. (coiih^- COS^) 



(5) 



the stress components can be expressed in elliptical coordinates in 
terms of a stress function Cf) as 



i 



j 2 sr 1 



sinhE COS h t 4-* 

7 a (sm^ -l-si-nh^) c)^ 



siril 7 cos~9 c) 4> 

l z (sin^ +sinh*|") c)rj 



( 6 ) 



I c) Z <£> _ sin77 cos*? c)^ 

3 2 drj 2 D Z (sin I ^ + sirih Z |') c)y 

+ smh cosh £ 

3* CSiri 1 ^ + Sibil f) SI - 



(7) 




7) z <i> sm7 cos^ ^ j> 
ctyDJ J 2 (sin^4si7ih^) c)l 



sin h ^ cos h £ c>4> 

J^sin^ + smh l §’) c>^ 



(8) 
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It is easy for one to verify that the expressions on page 
558 in [li] reduce to equations ( 6 ), (7) and (8) above. 

It might be noted here that the following assumptions have 
been made in the foregoing equations: 

i. Isotropic elastic material 

ii. Thin flat plate 

iii. Infinitesimal deformations 

The stress functions (c]a-j) assumed here are the same as 
those which Coker and Filon used in reference [H] 



C-n+i )§■ C-n-0i> 

Xl - e cos [(^- 07 ] + e cos 


[(^ + 1)7] 


(9a) 


X* - <2 cos^C'”- O' 7 ?] +e 00 s 




(9b) 


X3 — + e <n 


(9c) 


^ e <7 ’ +O? sm[('n-0^] + e. Cn 


(9d) 


X s = e cos ^*7 




(9e) 


Xg = e C 05 ^ 




(9f) 


-n^ 

X7 - e SK/t orry 




(9g) 
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Xa = e s'm Wfl 



(9h) 



X,= I 



(9i) 



X,o= *{ 



( 93 ) 



All the above functions satisfy the biharraonic equation for 
interger (or zero) values of or j. 

Equations (9a) through (9j) can be combined to give a general 
expression of the stress function 









00 1 


■ <n-n)$ 


(n- 


'hr' = 


°-°5 + °c 


+ „4 t ” 1 


e 


COS ( m - 1 ) 'y 4* 


e 




1 


r -Cn-0> 


cos (r»-l)rj 4. 


-<n- 


OS 






dn 


[ e 


e 


cos (n+i)i^ 


1 + 




1 


- 01+0$ 




fn - 1 


or 






«?n| 


e 


sin ( n~i)rj 4 


e 


sin (0 + 1)17 


1 + 






- -Cn+0? . 


-(r 


■w)| . 




+ 


fn | 


_ e 


Sin (n-i)T^ 4 e 


Sin (0+1)17 J 


+ 



O1 + O7] 



hn 



4 . gr) e COS nr} 4 . h n e cos nr} 4 . 



„ -it’ . 

4 - Ju e SmnT^ 4 -kn Q. Siring 



( 10 ) 
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When examinig the series with n=0 and n=l then for 

finite stresses and displacements it is evident that the 

coefficients d ,f ,h , j ,k ,f, must equal zero, 
o’o'o'-'o’o l 

The stress function given by equation (10) can be 
substituted in equations (6), (7), and (8) for the stresses. 

However in the present analysis , only loading conditions 
that maintain symmetry about both the x and y axes are 
considered and therefore the coefficients of the terms that 
produce nonsymmetric functions of y for different values of n 
were set equal to zero. 

Finally the stresses from (6), (7) and (8) have the form 



A r- r (nw)^ 1 

A2 1 % n f [ Cn ^ + '^ + On (rtfl) e Jcos(nHs)7 + 



4- 



(n+i)$ Cn-3)t -(n+Ot 

-4 Cm € - Cn Cm 3) e - A Dn € 



-C"- 3 ) 5 \ 

_pn(n-h 5 ) e COS (m Or] + 



+ 



Cms)!" 

4 Cn e - Cn (n-5) e e 



-"Dn (n-3) 



-(n+3)l” -r 

e 



Cos (n-0y 






Cn 



cmol" 



(n-l) e + Dn (n-0 e 



-Cn+0^ 



JcosCr7-s)y 



4 



- 15 - 



CO v V 

v — r r 

+ >n Gn(n-0« +Hn(n-i)e cos(n+ 2'>7 + 

nsj ,4 - ** 

cn-2')'^ -(n+a^'S 

+ ^-Gn(n-i)C -GnC ri +O e ~Hn(r>-i)e 



-(n-2)t 

- Hn (n+i)e ] a?s nrj 4 



+ [ Gn CTi+i)e n! + HnGn+ 0e nI ]cosCn-2^)y + 



+ 2 Ao s'mh 2 ? 



(ll) 



CO > 

. a \ r (n-0§ -(n-i )5 « 

A 3 6n-/_ T1 . [-Cn(n-3)e ~t>n(n-3)e JCO$ (ri43)y 4 



C A ‘ Tl = l, 3 ... 



4 [-4 Cn 



(041)? 



/=> 



* r. 



(n-J)? -<fn 4 0T 

(n-i)e ~ 4 Dn e 



+ Dn (n-|) e 



] cos (1141)7 4 



4 \A Cn 



e 



cn 43 )^ -(n-n^ 

4Cn(n+0C -+4Dn£ 4 



, -( 0 + 3 ) 1 " 

4Dn (n+i) e 



] C03Cn-l)7J7 



4 



Cn+i)^ 



[ V * 1 r 1 / O . 

~ Cn ( n + 3) & — Dn(n+3)<? 



-Cn + i)^ 



]co3(n- 



3)7 



4 



Z r 40 T -.nl* 

ni[-Gn(n-i)e -Hn(n-0e Jcos (71+2)7 4 

n=va,A .. v. 



- lH - 



(m2 






cn-*)! 






4 [Gn (n-OC- 4 Gn (n + i)6. 4Hn(n _ 0£- 4 



-(n-2.)t 

4 Hn (h+i)e J cosnT^ 4 
4[-6n(n4i) e 1 -Hn(n + 0 ~e ^]cos (tv 




- 2 Ao si n K 4 



( 12 ) 



oo 



43 T-* - 

"71" ' 

<- n = 1,3... 



n ^ [- Cn Cn-0r^Dn(n-0e ]siT7 (n43)l^-f 






4 [cn(n4i)e -T)n(n+0€. jsin (n4i')T^ 4 

r 0143)^ 1 _ 

4 [CnOvOe -PnOi-Oe JsmCn-Qr^ 






r £n-tO? - Cn+l)l’-i . . ] 

4 [-Cn Cn40e 4'Dr?Cn4i)e J s/n (n- s)rjj 4 

4 4 f [- Gr? (n-Q e 4 Hn (n-i) e j sin(n42)^ 4 



4 



Cn+z')? Cn-i)T -Cn 4 i)J 

Gn(n-Oe 4Gn(Vni)£ - Hn (n-i)e 

-(n- 2 )|' r 

-Hn (n4i) e jsin n-y 4 [- Gnfp4i)e. 4 



-nT 

4 Hn (n-u) C. 



SIT) 



(n-^)y A-Zho Sin 



(15) 
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However equations (ll), (12) and (15) are rearranged, in Appendix 1, 
in a more suitable form for translation into computer language. 

The normal and shear stress components in Cartesian coordinates 
can be computed from the equations ; 



G* - 6^ cosV + sincS sin.£°( (ih) 

Gy = sinV d cosV + T|^ sin^c< (15) 

T*y r 1 -G 7 ) si n 2 o< - cos ^ (16) 

co 5 o( — cos h ^ (l?) 

l^/cosh ^ -cos<?j> ' 

Sirio(— ^ sin_h_t C057 (18) 

[/ cosh. ^ -cos 2rj 



Equations (l4), ( 15 ) and (16) are derived as usual from 
equilibrium conditions on an element .(Fig, lb) 
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boundary conditions 



Equations (11), (12) and (13) satisfy the required conditions 
of symmetry and antisymmetry refered to in the previous section. 

In addition to this G| and T§p must vanish at . Page 193^ lj 

Therefore at each point on the elliptical hole we must satisfy 



Sj = o 




(19a) 


t ?T = ° 


& 

II 


(19b) 



For the case of one elliptical hole in an elastic finite 
flat rectangular plate loaded biaxially we must satisfy the following 
two conditions (2o) and (21 ), in addition to (19) 



T x\f r O (20a) 

<3>| - <oa ; . y- (20b) 



T. xy— O 



(21c) 



Kx ~ G i 



X 






(2! b) 
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Due to the symmetry about both the x and y axes it is 
enough to examine only one quadrant of the plate. 

Figure 2 illustrates all the cases under consideration. 

The numbers appearing in parentheses are the numbers of the 
equations of the boundary conditions applied for each particular 



case. 
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NUMERICAL METHOD 

Initially the number of unknown constants appearing in equations 
(|l) up to (13) was determined by an arbitrary truncation of the series. 

Later the behavior of the series was examined for several cases 
and it was realized that the series doesn*t always behave well. Similar 
convergence problems for functions involving elliptical coordinates were 
pointed out in reference [6] 

In this case, the reason is believed to be the very rapid convergence 
of the coefficients of Dn and Hn after the fifth or sixth term, resulting 
in a very rapid divergence of the coefficients of Cn and Gn for all of 
the cases that were examined. 

Finally twenty terms plus the constant term of the series were 
chosen. 

Using the Overdetermined Collocation method (Ref. 8 ) the following 
cases were examined: 

a. - Six fixed points on the inner boundary and another six fixed 

points on the outer boundary were selected to satisfy the 
boundary conditions. (Fig. 2d) 

b. ~ Decreasing by one the number of points on the ellipse with 

no change on the outer edges produced no difference in the 
numerical results obtained in case (a). 

c. - With six points on the ellipse and eight on the outer 

boundary, the shape of the curve for the tangential stress 
distribution around the elliptical hole was changed to a 
flatter one and it was realized that the boundary conditions 
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were not satisfied as well as before due to the least 
squares process involved in the numerical method, 
d.- The case of five points on the inner boundary and eight 
on the outer one was examined and it was found that 
the error introduced ■was smaller than in case (c), but 
not as small as in cases (a) and (b). (Verifying the 
satisfaction of the boundary conditions) 

It might be noted here that a smaller number of points than 
in case (a) was used by Christiansen [ 15 ] for a similar problem. 

For the case of axial loading parallel to the minor axis 
of the ellipse and for values of ^ < 0. 5 (wide plate) the 
expression approaches infinity and no results could be obtained 
(Appendix l) . 

The same effect happens for the case when plates with 3y>05 
are loaded parallel to the major axis of the elliptical hole. 

From the boundary conditions a set of M linear algebraic 
equations in IT unknowns is formed, where M > N . 

In matrix notation 



AX = B 



( 22 ) 



where A is the M by N matrix of the coefficients of 
the unknown constants X. 



By means of a least-squares method the normal equations 
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are obtained and the system in matrix notation has the form 



AAX = AB 



( 23 ) 



These normal equations are solved by elimination using 
largest pivotal divisor [Ref. 12 ] 

A computer program was developed and written in FORTRAN IV 
to execute all the required steps and calculate the tangential 
stress around the ellipse and the normal stresses along the 

axes of symmetry of the plate. 

Subroutines for matrix operation are taken from reference [ 12] 

All the necessary calculations for the problems under 
consideration are carried out by the same computer program by 
transfering control to the proper boundary conditions. 

The required input information for each plate to be 
examined, must be typed on one card and many plates can be 
analyzed the same time. 

Details for the input cards are given on comment cards in 
the program. 

Additional details about the problem formulation are given 



in Appendix 1 
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(a) 
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cc; 



T 

i 
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(d) 
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RESULTS 

Figures 5 up to 51 are typical plots of results obtained 
from the computer program. Whenever possible, experimental data 
are indicated. Most of the figures are self explanatory but additional 
information, if necessary, is given in the discussion of each figure. 

Figures 5 up to l4 show the stress distribution around the 
boundary of the elliptical hole in finite plates for several values 
of jj. and j ■ 

It should be noted that in cases of loading parallel to the 
minor axis, the tensile stresses predicted are 2 $ to 17$ higher 
than experiments reported in [ 9] while the much smaller compressive 
stresses are i(0$ to 60$ lower. 

When the external loading is parallel to the major axis of the 
ellipse the agreement with [ 9 ] lies between 5 $ and 25$. 

The above differences in compressive stresses, when loading 
parallel to the major axis are of minor importance from an engineering 
point of view because, these stresses are very small compared to the 
maximum tensile ones. 

Figures 15 up to 22 show the maximum tensile stresses (SCF and 
SCF') as function of j for several values of fu, . (Finite plate 
with one elliptical hole). 

In figure 25 several cases of biaxial loading are plotted. 

It was not possible to examine the same plate used in [10 1 because 
of n limitations discussed previously. Other plates were examined 



- 53 - 



and the 'definitions from [10] were used. 

The curves appear to have the same shape. 

Figures up to 51 show the Gy and Gy distribution around the 
boundary of the elliptical hole and the x axis respectively for typical 
uniaxial and biaxial loading. 



CONCLUSIONS 



The stress functions assumed have been shown to give results 
that are comparable with available experimental data and therefore 
some confidence exists in the theoretical predictions for cases 
where no experimental results exist. 

The stress concentration factors for plates perforated with 
elliptical holes under in-plane loading, appear to be much higher 
in comparison with similar plates perforated with circular holes. 

Certain disagreements have been noticed between measured 
and predicted stresses in the low level compressive stress range, 
while the predictions for the high level tensile stresses are 
somewhat higher than the measured ones but still in close agreement. 
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REC OMt'IENDAT IONS 

The following topics are recommended for further investigation: 

a. - Postulation of similar stress functions in terms of 

hyperbolic functions instead of positive exponentials, in 
association -with other functions decaying rapidly as 
increases in order to bo more suitable for computer use. 

Timpe in [15] suggests such stress functions. 

b. -i Study of the problem for cases of loading that produce 

antisymmetrical conditions (shear forces applied on the 
edges of the plate for example), will be rather easy to handle. 




<v 
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APPENDIX 1 

Equations (11), (12) and (15) may be rearranged in the form of (24), 
(25) and (26), respectively. Also equations (l4), (15) and (16) are 
rearranged in the form of (27), (28) and (29) respectively. 



%- 



43 \ Cn a 






('n-')l' 



Cn4i)T 



[0» + 0 c °s (*H3) 1 ^ 4- A cos (n-i)^ j 



4 



4 £ [-4 cosCn+0^ 4 (71-0 cos Cn - 3 M 4 

4 e [ "(n+3) c©S(V>+i>p 4 






4 e 



[ - Cn- 3) cos (h-/)op j j 



4 ) Datl 

n= 



-(n-i)t 



e [ Cn4i) cosCnt ?>)7 44 cos (n-i > 7 ] 4 



-Cn+O^. 

•4- e [ - 4 cos (ri4i)^ 4 (n-i) Cos (n-3)^j 4 

-fn-3Df 

4 € [ - (H4-5) Co sCo-^O^] 4 

_ (n+?oF 

+ <2 [ - (n-t>) cos (n-i)^ j | 



4 ) Gr] n 

0= *,4.. 



» (n-X)% 


n ! r 


<2 [ - Oi-tO cos ny ] 


4 e [ 



4 (n4Q cos (n-2)7^j4 
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Cn42)t 




"V + 



+ e [-Cn-Ocos 



r -(n-i)T - -n^ 

Hn n e [ -(n+o cos nij J + e [ (n-/) cosfn 45 )^ + 



-Cn +2) £ 

+ (n + i) cos Cri-2)7p] 4 € [ ~(n-0 cos >177 ] 



4 



4 A 0 si n zt 



(2k) 



^ — Cn-3)£ Cn-')^" 

S ±y -) Can \e [(n-OCtfSCn-tQ?] +e [ -(n-^ cos (n+3)y + 



? r a 

C n = i,).. 



Cn+Ol 



4 A cos C^-O^ ] 4 <2 [-4COS (n+i)^ - 



01+35^ 



- Cn-v 3) cos(n-3)r7 ] 4e [(n+OcDsCn-O?] J 



Z r -Cn-s) - 

t)n ii | e [(n-Qcos Cn+C>Y ] 4 [ -Cn-^)caS(n43)i^. 

n — I *• 



Cr>-Ot 






-(n+O^ 

-4 4 cos (n-O^ ] 4e [ - 4cosCr?+i)9p- 

-Cn-tS)^ - 

- (n+ 3 (> Cos (n- 3 )- 7 ? ] + e [Cn+i)coSCn-fV?] V 

Z 00 r <n-2)f nt 

Gait J e [(n+Q cos rm ] 4 £ - (n-i) cos Cn^z)^’ - 



Tlr3 ,4. .. 



^n+Z'tf 



- (n + 0 cos (n-2*)'?j?j +£ [ Cn-i)coSn7^] j 4 



CO 



4 /* Hn+l f e [( 77+0 cos ny] 4 C" [-(n-i) cos 07 + 2 ) 7 ? - 
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. . -.V 

- trt -*ij J 



— \ ■ ». * ** S « 

- (n+0 cos(n- 2 )r]\ +e [(n-i) cos rt'rfl | + 



Ao ^ — 2 ^ i n h 2^ ") 



(25) 



CO 



t A„n |e" A<"+O s| n(n+0$-e AnCrHS^V 

5 C* -n--!,!.... '• 



(n+0£ tr,+Vi , 1 

-e (n-fi) sin(n-3)r? 4- Q. Cn-i) sir?(f?-0^ j + 

xA- r -M ' -(n-o3> . 

4/ p n YH-e Cn+Osin(n40Y + e sin (n^Jy 4- 



3)$ 



4 € (n+OsinCn-a^-C. stn^n-0?7 ^ ^ 



-fn+3)^ 



Cn-z)X 



A 



4 ^ Gn Tl f e (>141) Sm -n 74 - £ (n+O sin(n-2^?4 






Cn+L)t ^ 

4 Cn-0 s'm Cn4-2)7^ 4 e (n-0 sir?7 ’2? j 



4 



^ ,fn-2)f -A . 

+/ Hn'D-j-e (n+O Si nn^fe (n4i)Sm (n-z)n? -f 

- 1 n +4)^ . - 

4 (n-0 Sin(n+4)r? 4 e (-n +0 si7?77r ? y 



U--M- • 



4 Ao 2 sin 2 



7 



( 26 ) 



Q0 



g A A _ y> Cn. smV Am 4 cos 1 ^ Bm - sin ^ ABm j 4 



mi,v.. 



1 
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4 y~t>n [sin % d P\xr\ 4 cos^ Ban - Sm 2 d Abzn] 4 

n- 

co 

Gn [ sine< A&n 4^osV"&3n — sin3?<f ABgn j 4 



<x> 

4- l-|n [ Sm^ A^n 4- Cos ^ p 

n - 2,4... 



— Si r\2U A^4n J 4 



4 A 0 [ - 2 sin 2 ^" cos 2 c< -2 sim 2 U sin ] (27) 



co 

A 3^ _ C n [ CO£*o{ Am 4 S 1 n V Bin 4 Sm 2^ AE>m ] 4 



00 



^ ^ B>n[ 'A^n 4Sin<2~Bzn -4Sin2^ ABzn J *4^ 



v 



4 - ^ Gn[ co'sd A$n 4 ‘Simx? "Ban -\-Sin 2 d A63n] 4- 



■n- 2 , 4 ... 



4 ^T~ Hn[ cos'd Cn 4 SinV^n 4 S in 2 d A&m] 4 



rtz 2,4... 



4 A 0 [ 2cosd sink 2^-2 sin d Sin 2 ^ 4 
42 sin 2d sin Zrj> 1 



( 28 ) 



oO 

Xy j 4 l\\ C — (Am - Bm') sin 2 d -ABin Cos 2n<|4 

C * 



c o 



+ > Cl 



n=i,5... 



j_ (Azn-'bm') sm 2d -AI5a n Cos 2d 



+ 



-6 o - 



0 ° 



+^> Can 


^ A jTi — Bd>r)^Si n 2o( — AB 50 CoS OoL J 






CP 


r _L ( A^n -Abse CDS 2d] 

* 2 


+)> H n j 



n^. 2,4... 



4 - 

4 



+ A 0 [(sin 2 ^ -v sinh 2^ sm3*f - 
- 2 sin cos 2c< 



( 29 ) 



In the above expressions for <Sx , Gy ; T*y : 

Ain (i = 1, 2, 3# ! 0 represent the coefficients of 

Cn, Dn, Gn, Hn, respectively, 
from the expression for <S? 

Bin (i = 1, 2, 3, *») represent the coefficients of 

Cn, Dn, Gn, Hn, respectively, 
from the expression for 6^ 

ABin (i = 1, 2, 3 , *t) represent the coefficients of 

Cn, Dn, Gn, Hn, respectively, 
from the expression for Z/y 
The value of is obtained by using the formula 



£= ±loq (.].+ b'/a- \ 

Z d / - b'fal } 



(30) 



-Gi- 



lt is obvious that as ^ -?• o (i.e. <x«W*) b’-va' and therefore 
— > oo . 

Analyzing the case of figure 2 d, for example, the use of 
equations ( 24 ) and (26) according to (19a) and (19b) generates a set 
of 12 equations, if six fixed points are selected on the boundary of the 
(quadrant) elli.pse 

Another set of six equations is generated by using (27) and (29) 

n 

according to (21), if three points are selected on x=3. 

Finally six more equations are generated by using ( 28 ) and (29) 
according to (20), if another three points are selected on y=3. 

It should be noted that the number of unknown constants was 
selected to be 21. 

The system obtained this way has 24 equations with 21 unknowns 
and it is solved by following the procedure described in section 
NUMERICAL METHOD. 

The matrix A T is formed and A and B are premultiplied by 
AT and from the solution of the system 



A T A X = AT B 



we obtain the values of the unknown constants. 

Summarizing the above procedure, we have chosen a number of 
points greater than the number of unknown constants and we required 
that the boundary conditions be satisfied as nearly as possible at 
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these points (Least square sense) Ref. [ l^t]. 

Because of the nature of the assumed series solution, the number 
of unknown constants is always odd when truncating the series at any 
value of n. Therefore the system of equations is always overdetermined 
at least by one equation. 
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